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文章 编号 :1005-3085(2009)06-0990-07 
具有 非 负 困难 度 的 符号 几何 规划 的 一 种 分 解 方法 * 


张 晓 梅 ， 王 燕 军 
(上 海 财经 大 学 应 用 数学 系 ， 上 海 200433) 


摘 要 : 本 文 针 对 具有 非 负 困难 度 的 符号 几何 规划 问题 提出 了 一 种 新 的 分 解 方法 。 该 方法 首先 利用 指数 变 

换 及 矩阵 理论 ， 将 原 问 题 等 价 地 转化 为 一 个 非 线性 程度 较 低 的 可 分 离 规划 ， 然 后 ， 将 所 得 等 价 问 
题 分 解 成 一 系列 易于 求解 的 子 问题 ， 并 且 当 困难 度 为 零 时 ， 文 中 给 出 了 求解 子 问题 精确 解 的 方 
法 。 最 后 ， 通 过 数值 实例 验证 了 新 方法 的 有 效 性 和 可 行 性 。 
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符号 几何 规划 作为 一 类 特殊 的 非 线性 规划 ， 已 广泛 应 用 于 产品 设计 、 任 务 管理 、 化 学 过 程 
设计 等 实际 问题 中 ， 为 最 优 设 计 提 供 了 有 力 的 工具 。 符 号 几何 规划 的 局 部 优化 方法 很 多 1 和 4， 
但 每 种 方法 都 有 各 自 的 缺点 ， 如 伪 对 偶 法 会 产生 严重 的 数值 困难 ， 压 缩 算法 虽然 优 于 对 偶 算 
法 ， 但 又 存在 收敛 速度 慢 的 缺点 。 针 对 这 种 情况 ， 本 文 在 文献 辐 的 基础 上 ， 提 出 了 符号 几何 
规划 的 一 种 分 解 方法 。 利 用 矩阵 分 块 理论 ， 将 原 问 题 等 价 的 转化 为 非 线 性 程度 很 低 的 可 分 离 规 
划 ， 在 求 出 此 问题 的 最 优 解 后 ， 利 用 原始 变换 就 可 得 到 原 问题 的 最 优 解 。 为 今后 求解 大 规模 的 
几何 规划 问题 提供 了 一 种 方法 。 

符号 几何 规划 的 标准 形式 为 目标 函数 为 单项 式 的 反 向 正 项 几何 规划 

min To 
(RPG): st. fml(£)<1, m=1,2,--,p 
fuz) 21, m —pctlp-ct2,..,M 
z0, 


其 中 
fm(7) = DI Cm ryen , 


Tm 是 函数 fis (m) 的 项 数 ，Cm: > 0，Ymtn 为 任意 给 定 的 实数 。fim(z) 为 正 项 式 。 


2 ”变换 与 分 散 
引 理 1 如果 对 规划 (RPG) 连续 作 两 次 变换 zn = e (n=0,1,…,N') 及 


EH oYmtnYn =logzmt, m=1,.……,M, 
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则 规划 (RPG) 有 解 的 充 要 条 件 是 存在 向 量 y 及 z， 使 下 面 规划 有 解 


min go 


St: ET ifs < 1, m = 1,2, "Dp 


XI Cauza zi m=p+1,p+2,..…,M (1) 
Ty = logz, 
z > 0， 
其 中 logz am (log zii: lgzum,- y ,logZmi,- log zmT t ^ ,log Zu, x ,log zu Ty)? > 
?110 Till sem YN’ 
^y120 ^y121 E ^'Y12N' 
工 一 
YMTm0O  "(IMTu1 ^^^ YMTMN’ 


证 明 (必要 性 ) EDAX (RPG) 有 解 ， 则 存在 z > 0， 满 足 


faz ) <1, m= 1,2,- ,p 
(2) 
2)21, m-ptlp-2,.-.,M. 


作 变换 ， 令 zn = et (n = 0,1,… NT), JHEAX (2) 中 ， 则 得 


XI m Call geminum = $m m Cmte? 3 < 1, Tn — 1,2, -P 
(3) 
Tn Ome genn E Ei "Cie N! o Ymtn n > 1, m=p+ 1, p +2,- ,M 


再 作 变 换 ， S OYminyn = log Zmt: 则 得 到 y 及 z 满足 (1) 中 不 等 式 。 又 因为 


yo e 
Zn=0711ngynr log z11 


yi 
Ty =T ` E : 一 : = log z 


EN oYM TmnYn log ZMTm 
Yn’ 


H zo = ew， 所 以 min zo 等 价 于 min yoo | 

(充分 性 ) 若 存 在 向 量 y 及 z， 使 规划 (1) 有 解 。 根 据 Ty = logz 知 zmt = eEn-oYmintn, F 
是 
PL Cmtzmt = Yr C, ae Pico, 
令 zn = 一 eyn， 则 结 直 论 成 立 。 

下 面 我 们 就 困难 度 d = 0 及 d > 0 两 种 情况 分 别 讨 论 规划 (1) 的 分 解 结 果 ， 其 中 困难 度 定义 
为 

d- XM T4—(N-1)$ S-—(N'41). 
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(I) a = 0 的 情形 
ENT A(N +1) x (N 十 1) 阶 和 矩阵， 不 妨 假 设 工 的 秩 为 (N’ +1) BEP É — I dE RE 
阵 。 于 是 有 y = 了 -11log z， 根 据 引 理 ， 规 划 (1) 就 等 价 于 如 下 具有 线性 约束 的 可 分 离 规划 


min eTT-!logz 


st. Eius <l, m=1,2, p (4 
XiCmz Rl m-optlpt2,.M 
z»0, 
其 中 e = (1,0,---,0)7 X (N + 1) ETUR. 
(I) 4d > 0 的 情形 
此 时 和 矩阵 工 的 行 数 大 于 列 数 ， 即 9 > N' + 1。 不 妨 假设 前 N’ 十 1 行 线性 无 关 。 设 
B 
r= ; 
Hop BJE(N' - 1) x (N' 十 1) 阶 非 奇异 和 矩阵， 相应 的 设 
log zg 
log z = : 
| log zp | 
由 Ty = logz 可 得 
Pd 
y= . 
D log zp 
根据 引 理 1， 规 划 (1) 等 价 于 如 下 规划 
min eT B-!logzg 
st. DB-llogzg —logzp 
bod oM UE < 1, m 一 12, m (5) 


Xi m m m-—pctlpct2,-.,M 


z » Q0. 


3 ” 子 问 题 的 求解 


(了 规划 (4) 的 求解 
Wer! = (u11, 212; ~t s UIT t ,UMI UM2)* ,UMTu )， 则 规划 (4) 的 最 优 值 可 通过 求 
解 下 列 M 个 子 问 题 获得 


min EP ug log zie min YT Upt log zpt 
(P1) : S.t. PR O21t Sa sm (Pp) : s.t. DN «1 


zu 0 Zm 0 
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min x in) p+ 1)t lOg z(p+ 1) min X EM um log zmt 
(Pp41) : s.t. per Cip+1)t?(p+1)t 21, ce , (Py) : s.t. SEN O HS 之 1 
Z(p+1)t > 0 ZMt >0 


在 此 基础 土 ， 我 们 分 三 种 情况 讨论 子 问 题 (B) 7 = 1,2,… ,p， 的 解 : 
(i) F uj = 0, t= 1,2,..….， Tj, 则 子 问题 (P) 的 最 优 值 为 Hj = (zh, 272，,"… o 2m) 是 
子 问题 (P;) 的 任何 可 行 解 ; 
(i) 车 we < 0, t= 1,2,… ,人 ， 其 中 至 少 有 一 个 不 等 式 严 格 成 立 。 
将 目标 函数 改写 为 
Xj ius log zjt = 2 log DE = —log (1 im. Ji 
由 于 II t1 2p 是 关于 Zjt 的 单调 增 函 数 ， H.zjc WEARI CjtZjt <1, 因此 满足 


T; zz 
3 Cj2jt zm] 


的 2, 为 最 优 解 。 
又 根据 文献 [6]， 
EL 
nu (sss £) —uji Ed E, ujt 
n2 2. 5t < Ut ZE Uje 
t-1^jt T; $a VO C T; n -ujt ? 
mala) ^ 0 mm) 
当 且 仅 当 
Cjzji..Cjazjiag _ Cjmjzjm; 
—Uuj —uj2 “UiT; 
时 等 号 成 立 。 故 子 问 题 (P;) 的 最 优 解 为 
Ujt 
zh = — i t= 1,2,- Tj, 
Cj ug 


最 优 值 为 
log (n. ($2) [or Ei), 


(iii) 若 至 少 存在 一 个 ua > 0， MS zito — 0, 有 Ujto log(zjto) > -œ TETN 

题 (P5) 无 解 。 

同样 ， 我 们 分 三 TRSRSOCHBUD pp M, BUR 

() 若 wjt = 0, t= 1,2,.… ,了 T), 则 子 问 题 (P ) 的 最 优 值 为 0， Hz; BR (P;) 的 任何 
可 行 解 ; 

(i) Buy 20,t—12,.,Tj, Hip —^oESCPERE HE. HTj-1iB BE 
12 — 二 ， 最 优 值 为 - Sb ET; > 2， 不 妨 设 wjs。 > 0。 È zito 一 0， 有 wjto log(zjto) 

一 一 co， "HET RIB (P) Xf. 

(i) EEDE Dut «0. WS > roo, 有 wjto log(zjto) > -oo， 于 是 子 问 
题 (P5) E. 

(ID 规划 (5) 的 求解 
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将 规划 (5) 中 的 约束 条 件 拆 成 两 部 分 ， 记 为 


min eTB-!logzg 
S.t. ET OmtZmt <1, m=1,2,...,p (6) 
ET” Cmtžmt > 1, m=p+1,p+2,...,M 


DB-! log zp = log zp. (7) 


设 eTB-1 = (ui,u5,:-- ,UN'+1)， 我 们 分 两 种 情况 讨论 规划 (5) 的 解 : 

à) 若 规 划 (6) 818. Bzg = ( 祝 , 经,… ia) zD = (zN+2 s) zp 为 (6) 中 约束 
的 任意 可 行 解 。 如 将 名 代入 (7) 中 ， 则 (7) 变 为 关于 log zp 的 线性 方程 组 ， 设 该 线性 方程 组 
RU zn. irpo (6) 可行， 则 原 规划 (5) 的 最 优 解 存在 ， 为 z* = (esp) HERE 
A DN tu; log zi; 

G) ”车 规划 (6) 无 解 ， 即 其 可 行 域 中 存在 一 列 {z}， 使 六 Twilog z*; > -oo。 若 同 
时 z* 也 满足 (7) 式 ， 则 原 规 划 (5) 无 解 。( 这 里 只 讨论 了 两 种 简单 的 情况 ， 其 它 情况 要 复杂 
很 多 ， 还 要 更 深入 的 研究 和 探讨 )。 


4 ”数值 实例 
例 1 考虑 以 下 问题 


min 224 z2 
sí. (Q0.3z1222 1, 


Z > 0. 
解 : 根据 引 理 1， 原 问题 可 等 价 的 转化 为 如 下 问题 
min —ilogzi-— 4 log z2 + log z3 
s.t. zi 29 € 1, 


0.323 2 1, 


z > 0. 


经 计算 d= 0， 于 是 上 述 最 优化 问题 等 价 的 分 解 成 两 个 子 问题 


min -i log z1 一 ilog Z2 min log z3 
(D: st z+z<l, (ID: st. 03221 
zı > 0, z2 > 0. 23 > 0. 


解 得 zi = za = l, za = 加 ， 于 是 原 问题 最 小 点 为 zf = cj = 1.8258， 最 小 值 夺 = 6.6665， 这 
与 文献 [1] 的 方法 得 到 的 结果 是 相同 的 。 
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例 2 考虑 以 下 问题 
min 2o pi? 
st. x5 gy <1, 


2 


—2,.—1 一 1 
Z1 73 十 Z273 2], 


zc. 
解 : 原 问题 可 等 价 地 转化 为 如 下 问题 
min — log zı — 0.8 log z2 — 0.4 log zs — 0.4 log z4 一 0.4log zs 
s.t. 2 <1, 
z2 +23 <1, 
Z4 + 25 21, 
z > 0. 


由 于 d = 0， 将 上 述 最 优化 问题 等 价 地 分 解 成 三 个 子 问题 


min -—logzi min —0.8log z2 — 0.4 log zs 
(D: st za<1, (I): st. z2+z3<1, 
2; > 0, 22250, z3>0, 


min 一 0.4log z4 — 0.4 log zs 
(IIT) : st. 2402521, 
z4>0, z >0. 


解 得 z1 一 1 Z2 = 2, Z3 = i Z4, 25 无 解 。 
例 3 考虑 以 下 问题 


min 120z1 + T2 
st. 到 zz <l, 
zi € 1, 
27120, 2zZ2 > 0. 
解 : 由 于 dg = 1 > 0， 故 原 问 题 等 价 的 转化 为 如 下 问题 
min —0.75 log zı — 0.25 log z2 — 0.25 log za 
s.t. 12025 +22 Xl, 
Bz <1, (8) 
z4 S1, 
2 > 0. 


0.25 log 2; — 0.25 log z2 — 0.25 log zs = log z4, (9) 
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由 (8) 解 得 
7 


p 45' 
YE zi, z2, 2 代入 约束 (9) 中 ， 得 zi=0.6344， 显 然 满足 约束 (8)。 所 以 2 zo, z3, za 是 原 问题 的 
最 优 解 。 代 回 变换 ， 得 原 问题 最 小 点 为 z+ = 0.6344, r4 = 25.1762， 最 小 值 zx = 100.7055. 


ZL = 0.0063, 29 = 0.25, 23 
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A Decomposition Method for Signomial Geometric Programming 
with Nonnegative Degree of Difficulty 


ZHANG Xiao-mei, WANG Yan-jun 


(Department of Applied Mathematics, Shanghai University of Finance 
and Economics, Shanghai 200433) 


Abstract: In this paper, a decomposition algorithm for signomial geometric programming with non- 
negative degree of difficulty is proposed. Through the exponential transformation and matrix theory, 
a problem equivalent to the original problem can be obtained which is a separable programming with 
lower nonlinear degree. This equivalent problem is decomposed into several subproblems, and the ex- 
plicit solution to each subproblem can be acquired when the degree of difficulty is zero. At last, some 
examples are given to verify the feasibility and efficiency of our algorithm. 

Keywords: signomial geometric programming; degree of difficulty; polynomial geometric program- 
ming 


